Interface currents in topological 
superconductor- ferromagnet heterostructures 



P M R Brydon^, Carsten Timm^ and Andreas P Schnyder^ 

^ Institut fiir Theoretische Physik, Technische Universitat Dresden, D-01062 
Dresden, Germany 

E-mail: brydonStheory . phy . tu-dresden . de 

^ Institut fiir Theoretische Physik, Technische Universitat Dresden, D-01062 
Dresden, Germany 

E-mail: carsten.tinmi8tu-dresden.de 

^ Max-Planck- Institut fiir Festkorperforschung, Heisenbcrgstrassc 1, D-70569 
Stuttgart, Germany 

E-mail: a . schnyderSf kf . mpg . de 

Abstract. We propose the existence of a substantial charge current parallel 
to the interface between a nonccntrosymmetric superconductor and a metallic 
ferromagnet. Our analysis focuses upon two complementary orbital-angular- 
momentum pairing states of the superconductor, exemplifying topologically 
nontrivial states which are gapped and gapless in the bulk, respectively. Utilizing a 
quasiclassical scattering theory, we derive an expression for the interface current 
in terms of Andreev reflection coefficients. Performing a systematic study of 
the current, we find stark qualitative differences between the gapped and gapless 
superconductors, which reflect the very different underlying topological properties. 
For the fully gapped superconductor, there is a sharp drop in the zero-temperature 
current as the system is tuned from a topologically nontrivial to a trivial phase. 
We explain this in terms of the sudden disappearance of the contribution to the 
current from the subgap edge states at the topological transition. The current 
in the gapless superconductor is characterized by a dramatic enhancement at low 
temperatures, and exhibits a singular dependence on the exchange- field strength 
in the ferromagnetic metal at zero temperature. This is caused by the energy 
shift of the strongly spin-polarized nondcgenerate zero-energy flat bands due to 
their coupling to the exchange field. We argue that the interface current provides 
a novel test of the topology of the superconductor, and discuss prospects for the 
experimental veriflcation of our predictions. 
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1. Introduction 

The discovery that gapped single-particle Hamiltonians can have a nontrivial topology, 
depending on their dimensionality and the presence of time-reversal and particle- 
hole symmetries [3 [51 [31 HI [5], has sparked a massive search for topological 
materials. A key motivation is to realize exotic Majorana-fermion states, which 
are guaranteed to appear at the edges of a fully gapped topological insulator or 
superconductor by the bulk-boundary correspondence [4] [5j, and which may have 
applications in quantum computation. Parallel to these developments, the concept 
of topological nontriviality has been generalized to gapless systems, such as nodal 
superconductors [i[71Hll|Tni[IIl[I21[I3j or Weyl semimetals [Hlllij. Bulk-boundary 
correspondences can also be developed in these cases, leading to the topologically 
protected appearance of nondegenerate zero-energy (or Majorana) arc lines or flat 
bands at certain surfaces. Much work has now been done on understanding the 
conditions under which these states can form, and a topological classification of stable 
Fermi surfaces of any dimension has recently been developed [131 116] . 

A promising class of materials in which to search for topological systems 
are noncentrosymmetric superconductors (NCS). The absence of bulk inversion 
symmetry in these compounds has two important consequences: it leads to a strong 
momentum-antisymmetric spin-orbit coupling and it permits the existence of mixed- 
parity pairing states with both singlet and triplet gaps present (17) . These exotic 
superconducting properties have inspired a strong research effort [TH], and many 
examples of unconventional superconductivity in NCS have been reported, e.g., 
CePtgSi pjj, CeRhSis [20], CelrSig [21], LiaPtaB [22], Y2C3 [23] and BiPd [24] ■ 
More recently, much attention has been focused on the possible nontrivial topology 
of NCS [S[i[Tni[IIl[Il[Tl[2Sl[Ml[171[lHl[Ml[3ni[3Ill32]- Specifically, the BCS 
Hamiltonian of an NCS belongs to class Dili of the Altland-Zirnbauer classification 
scheme. In two dimensions, gapped Dili systems may be topologically nontrivial 
and possess a nonzero Z2 topological number. An example of such a state is given 
by the NCS with Rashba spin-orbit coupling, s-wave form factor of the gap, and 
majority-triplet pairing [551 [5511271 [55] . As shown in figure [ija), the edge spectrum 
possesses helically dispersing subgap states with Majorana zero modes [5S1 [5^ [531 134) . 
as required by the bulk-boundary correspondence. In analogy to a quantum spin Hall 
insulator, the edge states carry a spin current [29l |34l [35] . Increasing the strength 
of the singlet pairing ultimately leads to a sign change of the negative-helicity gap, 
which marks the transition into a state with trivial topology [27) . The edge spectrum 
of this state does not display any subgap states. 

It is likely that many NCS are characterized by rather strong interactions, which 
may lead to higher orbital-angular-momentum pairing states )36) . e.g., the modulation 
of the gaps by a dxy-'Vfave form factor as shown in figure [Tj^b). Because of the presence 
of gap nodes, it is not possible to define a quantized global topological number for 
such an NCS. This NCS nevertheless also displays edge states with nontrivial topology. 
Every point in the (Im) edge Brillouin zone not lying on a projected gap node may 
be regarded as the edge of a fully gapped one-dimensional Hamiltonian in Altland- 
Zirnbauer class Am. This allows the definition of a momentum-dependent winding 
number W(^im), which can only change across projected gap nodes [9l[T0]. In the case 
of the (10) edge shown in figure [iTb), the winding number W(^iq^ evaluates to ±1 for 
ky between the projected edges of the spin-orbit-split Fermi surfaces. By the bulk- 
boundary correspondence for the one-dimensional class AIII Hamiltonian, it follows 
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Figure 1. The relationship of the bulk gap structure to the (10) edge states 
of an NCS with (a) s-wave and (b) d^y-wave form factors in a majority-triplet 
state. At the top of the figure we show the variation of the sign of the positive- 
and negative-helicity gaps, Ak,-|- and Ak,- respectively, about the spin-orbit split 
Fermi surfaces. Here red (blue) indicates positive (negative) sign, and we assume 
that the positive-helicity (negative-helicity) Fermi surface is circular with Fermi 
wavevector A:f,-|- {kp, — )- Below we show the spectrum at the (10) edge for singlet- 
triplet parameter q = 0.25, where black lines indicate edge states while the brown 
shaded region is the projection of the continuum onto the one-dimensional edge 
Brillouin zone. 



that there is a nondegenerate zero-energy flat band at these momenta [9l [30l |3T] . The 
dxy-'wave form factor is crucial here, as one may show that these states can only form 
when there is a sign difference between the gap on the forwards- and backwards-facing 
parts of the Fermi surface [TTJ|32]- As such, these states are present for both majority- 
triplet and majority-singlet pairing states. In contrast, the edge states at momenta ky 
lying between the projected edges of the positive-helicity Fermi surface depend upon 
the relative strength of the singlet and triplet gap: for a majority-triplet state there 
are topologically trivial dispersing states, whereas for majority-singlet pairing there 
are doubly degenerate zero-energy states with M^(io) — ±2. 

Much attention has been paid to the spin structure of the edge states of an 
NCS [23 [Ml [SS], as the surface spin current can be understood in terms of the 
polarization of the electronlike part of the edge state wavefunctions. We have 
recently shown [37) that the edge states of an NCS typically also have rather strong 
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total spin polarization, which consists of contributions from the electronhke and 
holehke components of the wavefunctions. Like the polarization of the electronlike 
wavefunction components [351 1311 i the total spin polarization is odd in the edge 
momentum as required by time-reversal symmetry, and depends on both the spin- 
orbit coupling and the relative strength of the singlet and triplet pairing. On the 
other hand, the total spin polarization can be rather different from the electronlike 
spin polarization, and it is the total polarization that couples to an external exchange 
field. In particular, we have demonstrated that the topologically protected zero-energy 
flat bands characteristic of the NCS with dxy-v/ave form factor have particularly strong 
total spin polarization. Coupling to an exchange field therefore gives opposite energy 
shifts to these flat bands on either side of the edge Brillouin zone, hence generating 
an imbalance between the integrated spectral density at these two momenta. This 
causes the appearance of a finite edge current, which depends rather strongly on the 
orientation of the exchange field, and shows a remarkable singular dependence on the 
exchange-field strength at zero temperature. In contrast, the absence of the flat bands 
for an NCS with s-wave form factor leads to a very weak edge current due to the 
interplay between the spin-orbit coupling and the spin polarization induced by the 
exchange field. 

The analysis in Ref. [37, was performed for an NCS strip in contact with a 
ferromagnetic insulator at one edge. The strip was described by a lattice model with 
an exchange field applied to the sites on one edge. Note that applying an exchange 
field to the entire NCS also leads to an energy shift of the edge states [131 EU IMI 
but additionally distorts the Fermi surfaces in a way which is inconsistent with a 
zero- momentum pairing state |39l 140] . In the present paper we consider an NCS in 
proximity contact with a bulk metallic ferromagnet (FM) , with each phase occupying 
a half-space and treated in the continuum limit. We note that such an NCS-FM 
heterostructure has been studied by other authors [4TJ |42], but they only address 
the proximity effect on the FM which does not concern us here. A fundamental 
difference between Ref. [37] and the present work is that the direct coupling between 
the edge states and the exchange field in the former is absent in the latter. While 
the edge states remain well-defined when an exchange field is applied to the edge 
layer [37], for the half-space continuum system tunneling into the bulk states of the 
FM gives them a finite lifetime and hence turns them into broadened resonances. One 
of our main goals is to understand how these differences affect the edge current. To 
accomplish this, we utilize a quasiclassical technique to express the interface current 
in terms of Andreev reflection processes. The quasiclassical technique also has the 
advantage of being able to treat realistically small superconducting gaps, whereas the 
exact-diagonalization approach used in Ref. |37j requires rather large gaps to avoid 
finite-size artifacts. We perform a systematic study of the current's dependence on 
the temperature, the pairing state of the NCS, and the exchange field in the FM. 
For the dxy-wave form factor, we find that the key features of the interface current 
are robust to the additional complications of a metallic FM. In contrast, the results 
for the s-wave form factor show that the broadening of the subgap interface states 
leads to qualitatively different behaviour of the current at low temperatures. We use 
the close relationship between the current and the interface local density of states to 
understand the origin of the current, and show how it reflects the topology of the 
NCS. 

Our paper is organized as follows. We commence in section [2] with the theoretical 
description of the system, including the construction of the scattering wavefunctions. 
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Figure 2. Schematic diagram of the NCS-FM heterostructure considered in this 
work. The NCS and FM occupy the x < and a; > half spaces, respectively, 
and are separated by an insulating barrier. The magnetization M of the FM is 
allowed to point in any direction. In general, we find that an interface current ly 
is present in the NCS as indicated by the black arrow. 



an ansatz for the Green's function in the NCS and the derivation of the current. 
In section [3j we present and discuss the results for the s-wave and d^y-wave gap form 
factors. This is followed in section |4] by a discussion of possible experiments. We 
summarize our work in section [S] 

2. Theory 

We study a planar junction between a bulk NCS and a metallic FM separated by a 
thin insulating barrier, see figure [2j The Hamiltonian is written as 

H = [ *^(r)i?(r)*(r) = / dV*^(r) ( foj""! ^j-^} ^ ] *(r) (1) 



At(r) -H^{r) 

where ^{r) is the Nambu spinor of field operators. We also introduce the notation of 
the hat and the caron (inverted hat) to denote 2x2 matrices in spin space and 4x4 
matrices in Nambu-spin space, respectively. The noninteracting Hamiltonian is 

do - A^s6(a:^)Hox ■ o" 



J?o(r)- EF + U6ix) 

I 2m 

+ ^ M-x) [(-iV) X e,] • CT + [(-iV) X e,] • ae{-x)} . (2) 

For simplicity, we assume that the effective mass m and Fermi energy Ep are the 
same in the NCS and FM. The insulating layer separating the NCS and FM is 
modeled as a (5-function barrier of strength U. In the NCS there is a Rashba 
spin-orbit coupling A, while the FM is characterized by an exchange field Hqx = 
I Hex I (cos 7? sin C + sinTysin^By + cos^e^); note that the spatially-dependent spin- 
orbit coupling in ([2| is written in symmetrized and thus Hermitian form. The spin 
degeneracy of the circular Fermi surfaces is therefore lifted in both the NCS and the 
FM: in the NCS we have positive-helicity (-I-) and negative-helicity (— ) Fermi surfaces 

with radii kF,± = kFi\/l + A2 =F A) where A = XkF/i2EF) and kp = yj2mEF /K^; 
similarly, in the FM the Fermi wavevectors for the majority (f) and minority spin {\.) 
Fermi surfaces are fc_F,^(^) = fcj^-^/l -I- (— )/icx7 where /lox = t^B\^cy\/ Ef- The minority 
Fermi surface disappears at /icx — 1, and the FM becomes a half- metal. 
The gap matrix is 

A(r) - Q{-x) /(V) i[i^ao + d • a] cjy (3) 
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where "ip = A(T)g and d — (— iVxez)A(T)(l— g)/fci? are the singlet and triplet pairing 
fields, respectively [T7]- The parameter q tunes the NCS from purely spin-triplet 
{q = 0) to purely spin-singlet (g = 1) pairing. The gap magnitude A(T) is assumed to 
have weak-coupling temperature dependence, with A(T = 0) = Aq. The form factor 
/(V) models different orbital-angular-momentum pairing states: /(V) = 1 for an NCS 
with (s-l-p)-wave pairing symmetry; and /(V) — —2dxdy/kp for a (c?2;j,-|-p)-wave 
pairing state. Writing the gap function in momentum space and adopting the helicity 
basis, there is only intra-helicity pairing with distinct gaps Aj. j- — /(«k)(-0i(; ± Idkl) 
on the spin-split Fermi surfaces. For the circular Fermi surfaces considered here the 
negative-helicity gap vanishes at q ~ qc with q^ — kp^- / [kp,- +kp) « 0.55, which may 
be regarded as the boundary between majority-triplet {q < qc) and majority-singlet 
[q > qc) pairing states. 



2.1. Scattering wavef unction 

The Bogoliubov-de Gennes equation for the quasiparticle states V'(r) with energy E 
is written in Nambu-spin space as 

H{y)iI){v) = EiI){y) . (4) 

We solve Q for the scattering states. As an example, the wavefunction for a i/-helicity 
electronlike quasiparticle with wavevector kj^ = {kv,xiky) incident upon the FM from 
the NCS is given by 



E 

v' = ± 



+ Q{x)e^^yy J2 {[0™e-^-- + d^>™e''=--]e(fcF,.-|fc,|) 

+ [c:>™e-''"- + d:,0™e-«"-] e(|fc,| - kp,^)} . (5) 

For a; < 0, the ansatz describes Andreev reflection of i/'-helicity holelike quasiparticles 
with wavevector kj^/ = {k^i^x,ky) and reflection probability amplitude a1^, and the 
normal reflection of z/'-helicity electronlike quasiparticles with wavevector 'k.^i = 
{—kv',x, ky) and probability amplitude . For cc > 0, the quasiparticle is transmitted 
as a spin-cr hole with wavevector kc^ = {—k„^x,ky) and amplitude c^^. or as a spin- 
(7 electron with wavevector k^- = {k^,x,ky) and amplitude df,^, respectively. Here 
we have made the standard assumption of E <^ Ep and thereby approximate the 
magnitude of the wavevectors for electrons and holes to be equal; relaxing this 
approximation is not expected to qualitatively alter our results. Note also that the 
momentum ky parallel to the interface is a good quantum number due to translational 
invariance along the y-axis. If \ky\ is larger than the Fermi momentum in a given 
scattering channel, only evanescent solutions in this channel are possible. These 
solutions are characterized by an inverse decay length into the bulk NCS (FM) of 
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The wave function ([5| is expressed in terms of the spinors for the NCS and FM. 
In the FM, the spinors for electrons (e) and holes {h) are 

e""'cos^, sin^, 0, 0^ (6a) 
-e-"'sin^,cos^,0, 0^ (66) 
0, 0, e'''cos^,sin0 (6c) 

0, 0, -e^^sin^, COS0 . (6d) 

In the NCS, the spinors for electronlike and holelike quasiparticles with momentum k 

and helicity u are given by 

1 ■ T 

</'?ki = — ^ (wk,!., -i^ie'^Mk,!., J^ie^^Sk.i-Wk,!., Sk.i.t'k.i/) (7a) 
V2 

1 ■ ■ T 

'/'/^k^i. ^ ("k,!., -i^ie'^Wk,!/, J^je'^Sk,iyWk,i., Sk.z^Mk.z^) (76) 
where k — kp^,y {cos 9, sin 9), Sk,,y = sgn(Ak,i/) and 



Uk.i. = \/ ^ (8a) 



f^k,. = _ A2 ^ . (8c) 
Evanescent solutions in the NCS are characterized by the spinors 

- ^ ^, 0, o)"^ (9a) 



e,ky ,i> 



V2 



^ ^0,0,1,-.%^^)^ (96) 



,NCS 



where k.v — ^ ky — klj^^. 

The reflection and transmission amplitudes in ([s]) are determined from the 
boundary conditions obeyed by the wave function at the NCS-FM interface. Firstly, 
we require that the wavefunction is continuous at the interface 

l}^ey{ky,X = Q~) = %l^ey{ky,X = Q'^). (10) 

To ensure the conservation of probability [13], the wavefunction must also obey the 
condition 

dxil^ey{ky,x)\^^^+ - dxil^eAky,x)\^^^- = cFr,® [zkpd-Q - iXkpay^ tpeuiky^x = 0)(11) 

where Z = Ukp/Ep is a dimensionless constant characterizing the strength of the 
insulating barrier. These conditions yield eight coupled equations for the probability 
amplitudes. 

The calculation for a z^-helicity holelike quasiparticle incident on the interface 
proceeds analogously. The reflection and transmission coefficients in this case are 
denoted by a superscript h, i.e., a'^^,, b'l^,, etc. 
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2.2. The Green's function 

Generalizing the method of Refs. [HI 133 |3S1 HZ], we obtain the following expression 
for the retarded Green's function GJij(jg(r, r'; i?) in the NCS as a 4 x 4 matrix in 
Nambu-spin space; 

''-^Gl,c^{ky,x,x'-E)e^^y^y-y^ 



2-K 



E 



ih?\/l + y kp COS 9^ 

"{ 



NCS i(k 



NCS MK.^+qu)x 



iNCS 



-e(fcF,,- |fcj,|)e(feF,.' - \ky\) 



■'h.k., 



{ 



-NCS \t 



Jky{y-y') 



(12) 



We present only the result for x' < x < 0, which is sufficient to obtain all quantities 
of interest; the Green's function for a; < a;' < has similar form. In constructing (12), 
we assume that |Ak,jy| = |Aj^ J = ^ky,u (vahd for the gaps considered here), and we 
hence introduce the notation — ilfc„,i/- For the Green's function it is necessary to 
include the energy-dependent corrections to the electron- and hole-wavevectors, i.e., 
ke,L',x ~ ki, x + and k^ u x ~ k^^x where 

m^ky,u 

qi' — / ■ (13) 

fiVl + A2 kp COS0^ 

Note that we neglect contributions from scattering into evanescent states in the Green's 
function ansatz, as enforced by the step functions in (12). 



2.3. Transverse interface current 

The currents are derived from the continuity equation. The charge-density operator 
/5(r) = — e X](T ^(t('')V'cr(i") obeys the Heisenberg equation of motion 



d 



1 



p{r),HN] + [p{r),Hp]) = -V • Je(r) - V • J,(r) (14) 



where Hn and Hp are the normal-state and pairing Hamiltonians, respectively. The 
commutators of these Hamiltonians with the charge-density operator correspond to 
the divergence of the so-called electronic and source current density operators, Je(r) 
and Js(r), respectively [351137]. The total current density is obtained by calculating 
the expectation value J(r) = (Je(i") + Js(r))- Only the electronic term contributes to 
the transverse current. 

After a lengthy calculation, the y-component of the charge current density in the 
NCS can be expressed in terms of the retarded Green's function as 



ieh 
Am 



27r 



d__d_ 
dy dy' 

ehk. 



-y5.)G5,cs(r,r';i?) 



2m 



e\ - 



^NCs(^J/' -^'j ^) 



(15) 



E—^iuir, 
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where Sx — Aia.g{axT—(yx)- Note that the first term in the brackets in (15) originates 
from the kinetic energy, while the second term is due to the spin-orbit couphng. The 
current density can equivalently be written as 

Jy{x) = - j ^ j dEi^^kyp{E,ky,x)+^-^p^{E,ky,x)^np{E) (16) 

where np{E) is the Fermi distribution function and 

p{E, ky,x) = -^Im (Tr {G^csl^i., E)}) (17) 

p-{E, ky,x) = -^Im (Tr {Gi,cs(fc>/, a;, x; E)Sx}) (18) 
are the energy- and momentum-resolved local density of states (LDOS) and a;-spin- 



resolved LDOS at distance x from the interface, respectively. Although (16 1 is of 
limited calculational value, it is useful for interpreting our results. In particular, we 
note that the interface current depends on the reconstructed electronic structure at 
the NCS-FM interface only through these quantities. 

From examination of the Green's function, we determine that four distinct 
scattering processes contribute to Jy{x): intra- and inter- helicity normal reflection, 
and intra- and inter-helicity Andreev reflection. All these contributions show 
exponential decay into the bulk NCS on the scale of the coherence length = 
hvp/iT^^o)- The interface current density due to normal reflection and inter-helicity 
Andreev reflection processes are, however, further modulated by rapidly oscillating 
factors with the length scales ^{2kp)~^ and ^(2Afci?)~^, respectively. As these length 
scales are much shorter than the total current contributed by these processes is 
negligible and we henceforth ignore them. We thus find the total interface current ly 
in the NCS to be 



.0 

ly ^ dx Jy{x) 

J —oo 



■ (19) 



In deriving (19), we find that the current in the ;/-helicity sector contributed by the 



spin-orbit coupling is exactly 2v\/kp ,j times that from the kinetic energy. 



A similar expression to ( 15 ) can be derived for the interface charge current density 
in the FM in terms of the FM Green's function G''^yi{r,Y';E). As in the NCS, the 
Green's function of the FM includes terms due to normal and Andreev reflection at 
the interface. Due to the absence of a pairing potential in the FM region, however, 
the latter processes only appear in the off-diagonal elements of the Green's function, 
and as such they give vanishing contribution to the trace in the expression for the 
interface current. Since the remaining contribution of normal reflection to the charge 
current density oscillates on the scale of the inverse Fermi momenta, we expect the 
integrated current in the FM to be negligible compared to that in the NCS. 



3. Results 



We find that an interface current appears in the NCS for a FM with magnetization 
components along the x- or z-axes, and reverses direction with the magnetization. In 
agreement with Ref. [37], there is no current for a y-polarized FM. In the following 
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we will consider the two cases of a magnetic moment pointing along tlic positive 2:- 
and z-axes. We present results only for spin-orbit coupling strength A — 0.2 and 
barrier strength U = SEp/kp- Note that the symbols in the plots of the currents 
are included to distinguish the curves and do not represent the only data points. In 
performing the Matsubara sum we utilize a frequency cutoff of 100 Aq. Our finite- 
temperature analysis is extended to zero temperature by replacing the Matsubara 



summation in ( 19 1 by an integral over the imaginary frequency axis. We express the 
current in units of eEp/2iTh, which is twice the magnitude of the edge current of a 
chiral p-wave superconductor [?51 05] ■ 

3.1. (s+p)-wave pairing state 

In figure |3] we plot the total interface current in the NCS for the (s-l-p)-wave pairing 
state as a function of temperature [panels (a) and (b)], singlet-triplet parameter 
q [panels (c) and (d)], and exchange-field strength /lex [panels (e) and (f)]. For 
weak to moderate exchange- field strengths the magnitude of the current in the 
(s-l-p)-wave NCS is always very small compared to the edge current of a chiral p- 
wave superconductor. At low temperatures it converges to a value which is weakly 
dependent on q in both the topologically trivial (g > qc) and nontrivial {q < qc) 
states, but has markedly larger magnitude in the latter. Just on the nontrivial side 
of the topological transition, we see that the current increases rather steeply at low 
temperature, note, e.g., the q — 0.5 curves in panels (a) and (b). Plotting the zero- 
temperature current as a function of q reveals a step discontinuity at q = qc, as seen in 
panels (c) and (d). Although the current in the topologically nontrivial NCS is roughly 
three times larger for a FM polarized along the z-axis than along the x-axis, in the 
topologically trivial state the two magnetization directions give comparable results. 
The current shows typical linear-response behaviour for small exchange- field strengths 
/lox- Approaching the half-metal limit hex = 1, we observe that the zero-temperature 
current appears to saturate for the x-polarized FM, but for a magnetization along the 
z-axis the current increases super-linearly with the exchange field. 

Insight into the origin of the current can be gained by examining the momentum- 
resolved quantity Iy{ky), defined as the even part of the ky integrand in ( 19 1. In figureH 



we plot the evolution of Iy{ky) with the singlet-triplet parameter q [panels (a), (bj] 
and the temperature [panels (c), (d)]. This reveals that the states at \ky\ < kF,+ 
are responsible for the sharp drop in the current at q = qc, and also for the large 
disparity between the currents for the x- and z-polarized FM in the topologically 
nontrivial case. In contrast, the momentum-resolved current at fcF,+ < j^yj < kp,- 
is almost independent of the singlet-triplet parameter for T = 0.1 Tc, except very 
close to q = qc, where the negative-helicity gap closes; indeed, in panels (a) and 
(b) the curves for q = 1, 0.75 and 0.25 are obscured by the q = Q result. At zero 
temperature all curves for q ^ qc are identical for kp^+ < \ky\ < kp_ (not shown). 
This can be understood by noting that due to the absence of positive-helicity states, 
varying q only scales the energy dependence of the Andreev reflection coefficients by 
the changing negative-helicity gap, which is irrelevant at zero temperature. On the 
other hand, the absence of the negative-helicity gap at q = qc completely suppresses 
Andreev reflection, and thus there is vanishing current. Note also that the current at 
kp^^ < \ky \ < kp _ is comparable for the x- and z-polarizations of the FM. 

The most interesting features of the current in the (s-l-p)-wave NCS are the 
zero-temperature step discontinuity at the topological transition {q = qc) and the 
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Figure 3. Interface currents in the {s+p)-wave NCS for a FM polarized along 
the X-axis (left column) and z-axis (right column). Current as a function of: 
(a), (b) temperature T at exchange-field strength hex = Ms|Hox|/i?F = 0-3; (c), 
(d) singlet-triplet parameter q for exchange-field strength feex = 0.3, where the 
vertical dashed line indicates the critical value qc\ (e), (f) exchange-field strength 
?tcx for singlet-triplet parameter q = 0.25. 



much larger current for the majority-triplet pairing state. The momentum-resolved 
results allow us to associate the jump a,t q — with a suppression of the current at 
\ky\ < kF,+ , which naturally suggests that the enhanced current in the topologically 
nontrivial state is due to the presence of subgap states. This is also consistent with the 
sharp low-temperature increase of the current near the closing of the negative-helicity 
gap at q — qc (see, e.g., the q — 0.5 curves in figures [3]( a) and (b)), as the subgap 
states can only contribute to the current at temperatures less than the minimum gap. 

To see how the edge states might carry a current, we first examine the analytically 
tractable limit of vanishing spin-orbit coupling and a purely triplet gap. In this 
case the Bogoliubov-de Gennes equation Q in the NCS splits up into two 2x2 
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Figure 4. Momentum-resolved interface current in the (s-l-p)-wave NCS for 
a FM polarized along the x-axis (left column) and z-&xis (right column) with 
exchange-field strength /lex = 0.3. (a), (b): The current for various values of the 
singlet-triplet parameter q at temperature T = O.lTc. (c), (d): The current for 
various temperatures T at singlet-triplet parameter q = 0.25. In all panels the 
vertical brown solid (green dashed) lines are the projections of the positive- and 
negative-helicity Fermi surfaces of the NCS (majority and minority spin Fermi 
surfaces of the FM), as shown in panel (a). 



systems with opposite chiral p-wave gaps, one for each Cooper pair z-spin orientation 
s = ±1. At a vacuum edge of this so-called helical superconductor one finds subgap 
states with dispersion Es=±{ky) = sAoky/kp [HI ESI [M] ■ When placed next to a 
z-polarized metallic FM, the spin-dependent reflectivity TZs{ky) leads to an unequal 
broadening of these states. In particular, by calculating the poles of the Andreev 
reflection coefficients, we obtain the energy of the subgap states 



Es{ky) = Ao sin |^arcsin + i\n ^Jns{ky)j . (20) 

The left- and right-moving subgap states therefore have finite but different lifetimes 
when TZ+{ky) ^ TZ-.{ky) < 1, and thus appear as resonances of different width in 
the energy- and momentum-resolved LDOS p{E, ky, x) in the NCS. On the other 
hand, there is no asymmetry in the spin-degenerate continuum states, i.e., we find 
that p{E,ky,x) 7^ p{E,—ky,x) only for energies \E\ < Aq. Furthermore, the 
integrals of the LDOS at ky and —ky over the negative-energy states are unequal 

/"ao ^y^ ^) 7^ /-Ao Pi-^^ ~^yi ^^"^ so from ( 16 ) we deduce that there is 

a finite charge current density in the NCS which is carried entirely by the subgap states. 
We note that the spin-dependent broadening is not connected to the total polarization 
of the states, but rather is controlled by the polarization of the electronlike part of the 
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Figure 5. Zero-temperature energy- and momentum-resolved LDOS in the 
{s+p)-wave NCS at the interface with the FM. We present typical results for 
both X- and 2-FM polarizations with /lex = 0.3, and also both the majority- 
triplet and majority-singlet pairing states. The left column shows the LDOS, 
while the right gives the momentum-antisymmetrized LDOS Ap{E,ky,x) = 

{p{E,ky,x) — p{E,—ky,x))/2, both normalized by po = ml-KY?kp\/ 1 -f A^. We 
assume an intrinsic broadening of 5 X 10~'*Ao. 
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wavefunction. It follows that for the helical superconductor the broadening of the two 
edge states is the same when the FM is polarized along the x- or y-axis, and there is 
hence vanishing interface current. 

It is not possible to rigorously make the above argument in the presence of spin- 
orbit coupling or of a singlet gap. We nevertheless expect that the spin-dependent 
broadening of the edge states should be robust to these complications, as the spectrum 
evolves continuously as they are switched on. These spin-mixing terms may also 
lead to a similar broadening effect for an FM polarized along the x-axis, as they 
lift the degeneracy of the positive- and negative-helicity states [53]. To test this, 
we plot in figure [5] the energy- and momentum-resolved LDOS at the interface for 
systems representing both topologically nontrivial and trivial states and also both 
magnetization directions. An asymmetry in p(E,ky,x = 0~) is readily visible only 
for the case of a majority-triplet pairing state and a z-polarized FM [panel (e)]: the 
right-moving subgap state is clearly less broadened than the left-moving state. Subtle 
changes in the LDOS are nevertheless present for the other cases, which can be revealed 
by plotting the momentum-antisymmetrized quantity 

ApiE,ky,x^O-) = ^ [p{E,ky,x^Q-)- p{E,-ky,x = 0-)] (21) 

in the right column of figure [5j The different broadening of the two subgap states for 
the q = 0.25 NCS in contact with a z-polarized FM [panel (f)] now becomes much 
clearer. On the other hand, the results for the q — 0.25 NCS in contact with an 
x-polarized FM [panel (b)] reveal that the ky > and ky < states have been slightly 
shifted to higher and lower energies, respectively. This is consistent with the coupling 
of their spin polarization to the exchange field |37j : a similar but smaller energy 
shift for the z-polarized FM is masked in figure [5f^f) by a much greater broadening. 
Closer inspection of the subgap states in figure [sj^b) nevertheless reveals that the right- 
moving state is less broadened than the left-moving state, although the difference 
in linewidths is much smaller than for the z-polarized FM. For the majority-singlet 
cases, the antisymmetrized LDOS reveals that the changes in the LDOS are small and 
restricted to the gap edges. 

This analysis shows that the unequal broadening of the right- and left-moving 
edge states is a plausible explanation for the enhanced current in the topologically 
nontrivial state, and also why there is a much larger current for the z-polarized 
FM. Furthermore, it explains why a sharp jump in the current across the topological 
transition was not anticipated by the exact-diagonalization calculations presented in 
Ref . [37] . The broadening of the edge states is dependent on an imperfect refiectivity, 
so that spectral weight can leak into the FM. Clearly this is only possible for a metallic 
FM; furthermore, in our ballistic-limit calculation we require that the FM be of width 
much larger than the NCS's coherence length, so that electron and hole pairs decohere 
before being refiected back at the opposite side of the FM towards the superconductor. 
In contrast, in Ref. [37] an insulating FM was considered, which obviously preserves 
the perfect edge reflectivity, and indeed no broadening of the edge states was found. 
The resulting edge current therefore has completely different characteristics to that 
found here, as it is due entirely to the induced x-spin polarization via the spin-orbit 
coupling. 
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Figure 6. Interface currents in the {dxy+p)-^sve NCS for a FM polarized along 
the X-axis (left column) and z-axis (right column). Current as a function of: (a), 
(b) temperature T for exchange-field strength h^,^ = 0.3, with insets showing low- 
temperature T < 0.1 Tc currents in more detail; (c), (d) singlet-triplet parameter q 
for exchange-field strength hox = 0.3, where the vertical dashed line indicates the 
critical value q^; (e), (f) exchange-field strength hex for singlet-triplet parameter 
q = 0.25. 



3.2. (dxy+p) -wave pairing state 

As shown in figure|6j the currents in the {dxy-\-p)-'weive NCS are dramatically different 
from those in the (s-l-p)-wave state. We first note that the sign of the currents is 
reversed between the two cases, originating from the additional tt phase shift acquired 
by Andreev-refiected quasiparticles due to the d^y form factor. In further contrast to 
the (s-l-p)-wave case, the currents due to an x- and z-polarized FM are quantitatively 
very similar. More remarkable is the temperature-dependence of the current [panels 
(a) and (b)], which is characterized by a sharp increase at low temperatures for all 
q Qc] for q = qc, in contrast, the current grows only slightly below T = 0.1 Tc, see 
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the insets of panels (a) and (b). As seen in panels (c) and (d), the magnitude of the 
zero-temperature current for q ^ qc far exceeds that in the (s+p)-wave case. The 
zero-temperature current is discontinuous at qc- exactly at this point it takes a value 
roughly an order of magnitude smaller than at g = ± 0"*^ . There also appears to be 
a small jump in the current between q on either side of q^ accompanied by a change 
of the slope. The difference between the current in the majority-triplet and majority- 
singlet regimes is nevertheless much smaller than for the (s-|-j5)-wave NCS. At low 
nonzero temperatures, however, the current is sharply suppressed as one approaches 
q = qc, indicating that the low temperature enhancement becomes increasingly sharp 
near the negative-helicity gap closing, e.g., see the q = 0.5 curve in the insets of panels 
(a) and (b). Indeed, the absence of the low-temperature current enhancement for 
q — qc implies that it crucially involves the negative-helicity states. 

The dependence of the current on the exchange field [figures ^e) and (f)] 
shows a remarkable deviation from linear-response behaviour at low temperatures: 
for T < 0.01 Tc the current grows very rapidly with the exchange-field strength, 
before going through a maximum and slowly decreasing. At zero temperature an 
infinitesimally small exchange field in the FM is sufficient to generate a large current 
in the NCS. We note that there appears to be a qualitative change in the dependence 
of the zero-temperature current on the exchange field at /igx ~ 0.5, characterized by 
a change of slope of the low-temperature currents. 

The momentum-resolved current (figure [?]) shows that the low-temperature 
enhancement of the current is due entirely to states at kp^^ < \ky\ < kp^-, consistent 
with the critical role of the negative-helicity gap deduced above. Although there is 
considerable variation of this current as a function oi q for T = 0.1 [figures [7]ja) 
and (b)], at zero temperature there is no variation with q for q ^ qc (not shown), 
for the same reasons as in the (s-l-p)-wave case. Thus, the change in the total zero- 
temperature current across the triplet-singlet boundary [figures ^c) and (d)] is due 
only to states at \ky\ < kp^+. The temperature-dependence of the kp^+ < \ky\ < kp- 
current is astonishing: comparing the T = 0.01 Tc and T = curves in figures [Tf^c) 
and (d), we observe that while the current at < kp^^ has almost saturated to 
its zero-temperature value by T = 0.01 Tc, for kp^^ < \ky\ < kp- the current more 
than doubles as the temperature is lowered. For q = 0.25, the the zero-temperature 
current at fc^.t < \ky \ < kp- accounts for more than 40% of the total. Interestingly, 
in this momentum range the current displays a linear dependence upon ky, and is 
almost independent of the exchange-field strength [figures ^e) and (f)]. The current 
at kp,^^ < |fcj,| < kp,- is therefore clearly somewhat special, and its disappearance when 
fcF,t = kp,- at the exchange-field strength /lox ~ 0.49 closely matches the qualitative 
change in the exchange-field dependence seen in figures |6][e) and (f). In closing, we 
note that although the momentum-resolved currents at \ky \ > kp are very similar for 
the two polarization orientations, clear differences are seen at smaller momenta. 

The remarkable temperature dependence of the current in the ((ia;y-|-p)-wave NCS 
is intimately connected to the coupling of the exchange field to the topological edge 
states at < \ky\ < kp-. In Ref. [37] it was shown that at a vacuum edge these 
states possess strong x- and z-spin polarization, with equal magnitude but opposite 
sign at ky and —ky, as required by time- reversal symmetry. Applying an exchange 
field to the edge of the NCS therefore shifts these states in opposite directions, one 
above and the other below the Fermi energy. From (16 1 we see that the first term 
in the zero-temperature momentum-resolved current is proportional to the difference 
in the number of states below the Fermi energy at ky and —ky, while the second is 



Interface currents in topological superconductor-Jerromagnet hetero structures 



17 



X 




0.25 0.5 



1.25 0.25 
klk^ 

y F 



1.25 




Figure 7. Momentum-resolved interface current in the {dxy+p)-viSMe NCS for a 
FM polarized along the x-axis (left column) and 2-axis (right column). We show 
the momentum-resolved current: (a), (b) for various values of the singlet-triplet 
parameter q at temperature T = 0.1 Tc and exchange field /lox = 0.3; (c), (d) 
for various temperatures T at singlet-triplet parameter q = 0.25 and exchange 
field hex = 0.3; (e), (f) for various exchange fields hex at zero temperature and 
singlet-triplet parameter q = 0.25. In all panels the vertical brown solid lines are 
the projections of the positive- and negative-helicity Fermi surfaces of the NCS. 
In (a)-(d) the projected majority and minority Fermi surfaces of the FM are given 
by the vertical green dashed lines, while in (e) and (f) they are indicated by the 
colored arrows. 
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Figure 8. Zero-temperature energy- and momentum-resolved LDOS in the 
{dxy+p)-vSL^e NCS at the interface with the FM. We present typical results for 
both X- and z-FM polarizations with /lox = 0.3, and also both the majority-triplet 
and majority-singlet pairing states. The left column shows the LDOS for momenta 
\ky \ < lAkp and energies |_E| < 1.2Ao, while the right column shows a low-energy 
detail for ky > 0. In the latter, the light green and dark magenta lines indicate the 
location of the peak in the subgap LDOS and its FWHM, respectively. In panel 
(b) we indicate the projection of the majority- and minority-spin Fermi surfaces 

in the FM. In all panels the LDOS is normalized by pQ = m/irh^kp^/ 1 + A^, and 
we assume an intrinsic broadening 5 X 10^^ Aq. 
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proportional to the sum of the x-spin polarization at these momenta. The opposite 
energy shifts of the oppositely polarized flat-band states at ky and —ky thus can 
contribute a large current. This implies the breakdown of linear-response behaviour 
at zero temperature, as then even an infinitesimal energy shift of these states causes 
a discrete change in the number difference and the x-spin polarization. 

For the device studied here, this argument can be directly applied for kp^^ < 
\ky \ < kp,- where only evanescent solutions are possible in the FM and so the subgap 
states in the NCS have zero broadening. As seen in the LDOS plots in figure [8j 
these states acquire an energy shift < O.OIAq from the coupling to the exchange field. 
From (16), we find that the zero-temperature momentum- resolved current for these 
states has the form 

^^^''^^^^^fct?^' ^^^^ 
in the hex ^ limit. The induced polarization at finite values of the exchange 
field can slightly reduce this current, as seen in figures ^e) and (f). Although the 
momentum- resolved current ( 22 ) vanishes with the positive- helicity Fermi surface as 
the spin-orbit coupling strength diverges, the momentum-integrated zero-temperature 
current nevertheless saturates at eEp/Anh in this limit. Much larger values of the 
low-temperature current than presented here are thus theoretically possible. Because 
of the small correction to the energy of the flat-band states, their contribution to the 
current only develops at extremely low temperatures T < 0.01 Tc, which is indeed 
observed in figures [tJ^c) and (d) . We note the apparent paradox that in spite of the 
strong spin-polarization of these states at \ky \ = kp,^ [37], the LDOS plots show that 
the energy shift vanishes as \ky \ — >■ kp-, indicating that the states are "anchored" to 
the node of the negative-helicity gap. This is consistent with the observation that the 
current saturates at ever lower temperatures as one approaches ±kp^- (not shown). 

The flat-band states at kp^^ < \ky\ < kp^^ also acquire an energy shift due to the 
coupling to the exchange field, but the presence of the open scattering channel in the 
ferromagnet also gives them a finite lifetime. This could not be anticipated by the 
analysis in Ref. [37], as discussed in section 3.1 above. In the right column of figure [sj 
we mark the maximum in the subgap LDOS by the green lines; the broadening is 
quantified by the full width at half maximum (FWHM) curves shown in magenta. 
The increasing broadening as \ky\ — ^ fc_F,+ results in a rapid suppression of the zero- 
temperature momentum-resolved current, as the imbalance between the integrated 
weight at ~ky and ky is reduced. On the other hand, the broadening allows the 
current to saturate at temperatures well above the energy of the subgap maximum. 
Remarkably, it still seems possible to think of the current in terms of the energy shift 
of the subgap state, even when this is much smaller than the broadening. For example, 
the zero-crossing of the subgap maximum at \ky \ ss kp,^ for the x-polarized FM is nicely 
correlated with a sign change in the momentum-resolved current [see figures [7]ja) and 
(e)], as expected from treating the subgap maximum as an unbroadened state. We 
note that there is no zero crossing of the subgap maximum for the z-polarized FM, 
and also no sign change in the current. 

In closing, we note that the LDOS also shows significant structure for states lying 
within the projected positive-helicity Fermi surface {\ky\ < kp^^), i.e., broadened 
dispersing states for majority-triplet pairing, and zero-energy states for majority- 
singlet pairing. Although in the latter case there is a slight shift in the location 
of the subgap maximum, the large broadening washes out any contribution to the 
current, except perhaps for \ky\ very close to A;f,+ . 
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4. Experimental prospects 

The interface currents in the NCS-FM junction discussed above contain clear 
signatures of the topological state of the NCS. For the (s+p)-wave NCS this is the 
sharp jump in the current at the topological phase transition, while the presence 
of the topologically protected flat bands in the ((ia;y+p)-wave NCS directly leads 
to the sharp increase of the current at low temperatures and a nonperturbative 
dependence on the exchange-field strength. The detection of any of these effects 
would therefore be strong evidence for a nontrivial topology of the NCS. The current 
characteristics of the {dxy+p)-^s^re NCS are of particular interest as they only arise 
from nondegenerate flat bands. In contrast, previous proposals to evidence the flat 
bands by tunneling conductance measurements cannot easily distinguish between 
nondegenerate and doubly degenerate states [HI [30l [SU [32] . 

The experimental verification of our predictions is nevertheless likely to be 
challenging, and must overcome a number of obstacles. Foremost is the Meissner 



effect, as screening in the NCS will exactly compensate the interface current (19) for 
the half-space geometry considered here. For an NCS of finite width W , however, 
this problem can be avoided by exploiting the different length scales of the interface 
and screening current densities, the coherence length and the penetration depth Al, 
respectively. The best NCS topological superconductor candidates are extreme type-II 
superconductors, e.g., CePtaSi has Ginzburg-Landau parameter k = Xl/£,o w 140 [TB]. 
For such a material it is possible to choose the sample width such that S^q <^ W <^ Xi^, 
which implies that while our calculation for the edge currents remains valid, screening 
currents are very small. The penetration of the FM's magnetic field into the bulk NCS 
must then be considered: although the NCS pairing state is robust to weak fields along 
the z-axis, which are more easily kept outside of the NCS in any case, a field along the 
X-axis can destabilize the NCS towards a phase where the Cooper pairs acquire a finite 
momentum (33 HO]- Because the consequences of this for the edge states is unknown, 
a z-polarized FM is therefore a more favorable choice for experimental study. 

The construction of an NCS-FM heterostructure device presents further 
difficulties. Firstly, some degree of surface roughness is unavoidable, which will 
lead to additional broadening of the interface states [SO]. As long as this does 
not introduce further energy shifts, however, we expect that the interface current 
should persist for weak disorder. Indeed, we have seen above that even the strongly 
broadened edge states of the {dxy+p)-wa.ve NCS contribute a significant current. 
A more serious problem is the choice of material for the NCS part of the device. 
Although there are many examples of bulk NCS, little work has been done on 
incorporating them into heterostructures. An alternative approach is to instead 
engineer the NCS in the heterostructure, say by coating a superconducting substrate 
with a thin normal layer of a material with strong spin-orbit coupling, so that the 
former induces a superconducting gap in the latter. This is in the same spirit as the 
well-known proposal to artificially create a topological superconductor in a quantum 
wire [53 155], which has been reported in recent experiments [53]. All suggestions 
along these lines [5TJ [52 [5U [551 [5^1 [571 [SS] involve the modification of the standard 
proximity-induced superconductivity by spin-orbit coupling of the same form present 
in a bulk NCS. Using unconventional c?-wave cuprate |57j or s±-wave pnictide |58j 
superconductors holds particular promise, as it might then be possible to artificially 
create the most interesting cases of {dxy+p)-'waye or topologically nontrivial (s+p)- 
wave NCS, respectively. No matter how the NCS-FM heterostructure is constructed. 
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however, there will be some variation ol the superconducting gaps close to the interface 
due to the pair-breaking effect of the tunneling barrier and the FM [SS] ■ Although our 
calculation does not account for this, the current predicted above should be robust as 
it ultimately arises from the spin structure of the bulk condensate. 

5. Summary 

In this paper we have used a quasiclassical method to study the properties of the 
charge current that appears at the interface between an NCS and a metallic FM in 
a two-dimensional junction, where each phase is assumed to occupy a half space. 
We have considered two complementary models of the NCS: a gapped (s-fp)-wave 
pairing state, and a gapless (c?2,j,+p)-wave system. Due to the contrasting topological 
structure of the two models, we find completely different dependences of the interface 
current on the temperature, the exchange-field strength /icx and the singlet-triplet 
parameter q. In both cases we find signatures of the topology in the interface transport. 
For the (s-l-p)-wave NCS, the topological transition from the nontrivial to the trivial 
state is signaled by a discontinuous drop in the zero-temperature current, due to the 
disappearance of the contribution from the subgap states. In the (d2,j,+p)-wave NCS, 
there is an enormous enhancement of the current as the temperature approaches zero, 
and the dependence on the exchange-field strength becomes singular. This anomalous 
behaviour originates from the energy shifts of the spin-polarized flat bands due to 
the coupling to the exchange fleld in the FM. While the results for the {dxy+p)-^ci.'ve 
NCS were anticipated by the analysis of an exchange field applied directly to the 
edge [37], the current in the (s-l-p)-wave NCS is qualitatively different due to the 
broadening of the subgap states by tunneling into the metallic FM. We thus find 
that the mechanism based on fiat bands, relevant for the {dxy+p)-^scve case, is rather 
robust and independent of the detailed nature of the FM. Hence, one can speculate 
that similar current characteristics might also be realized in other systems possessing 
topologically protected nondegenerate flat bands. 
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